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Abstract. We present a theoretical investigation on pe- 
riods and amplitudes of RR Lyrae pulsators by adopting 
stellar parameters which cover the range of theoretical 
evolutionary expectations. Extensive grids of nonlinear, 
nonlocal and time-dependent convective RR Lyrae enve- 

^^O lope models have been computed to investigate the pul- 
sational behavior in both fundamental and first overtone 
fj modes at selected luminosity levels and over an effective 
^ temperature range which covers the whole instability re- 

^"5 gion. In order to avoid spurious evaluations of modal sta- 

lf) bility and pulsation amplitudes, the coupling between pul- 
sation and convection was followed through a direct time 

( «— ^ integration of the leading equations until radial motions 
approached their limiting amplitude. 

Blue and red boundaries for pulsational instability into 
\^ the HR diagram are presented for three different mass 
(~*> values M=0.75, 0.65 and 0.58 M & , together with an atlas 
of full amplitude theoretical light curves for both funda- 
CTN mental and first overtone pulsators and for two different 
assumptions of stellar masses: M=0.75 and 0.65 M & 1 . 
O jThe dependence of periods on stellar parameters is dis- 
l cussed and new analytical relations connecting the period 
O to the masses, luminosities and effective temperatures are 
^— > provided. 

^ We show that theoretical expectations concerning mini- 
mum fundamental periods are in good agreement with the 
observational evidence of a dichotomic period distribution 
between different Oosterhoff type clusters. A rather good 
correlation has been found between the pulsational ampli- 
tude of fundamental pulsators and the effective tempera- 
ture, rather independently of stellar mass and luminosity. 
Theoretical periods have been combined with theoretical 
amplitudes in order to predict the location of the pulsators 

Send offprint requests to: G. Bono, Osservatorio Astronomico 
di Trieste, Via G.B. Tiepolo 11, 34131 Trieste, Italy 
1 The Tables 5 to 10 which summerize the pulsational am- 
plitudes of the models computed in this investigation are only 
available in electronic form at the CDS via anonymous ftp 
130.79.128.5 



in the Bailey amplitude-period diagram. Comparison with 
observational data brings to light what we regard as a clear 
indication that the OR region, i. e. the region where both 
fundamental and first overtone show a stable limit cycle, 
is populated by fundamental or first overtone pulsators in 
Oosterhoff I and Oosterhoff II clusters respectively. Some 
evident mismatches between theory and observation have 
also been found, and they are presented and discussed. 

Key words: stars: variables: other (RR Lyrae) - stars: 
evolution - stars: horizontal-branch - globular clusters: 
individual (M3, M5, M15) 



1. Introduction 

The pulsational properties of RR Lyrae variables of galac- 
tic globular clusters (GGCs) provide an unique oppor- 
tunity for testing the structural parameters of pulsating 
stars and the prescriptions of stellar evolutionary theories. 
According to such an evidence, in the last decades the 
well known van Albada & Baker (1971) relation connect- 
ing the pulsation periods to stellar masses, luminosities 
and effective temperatures has been at the basis of sev- 
eral investigations on the evolutionary status of HB stars, 
raising in particular a well known debate about the actual 
luminosity of RR Lyrae pulsators. 

A relevant point to remember is that pulsational periods 
are unaffected either by interstellar reddening or by the 
distance modulus and thus provide firm and accurate ob- 
servational data which unambiguously constrain the evo- 
lutionary parameters of the evolving variables. However, 
the pulsational phenomenon provides a further relevant 
observable, namely the amplitude of the luminosity varia- 
tion. According to such plain evidence, as early as the be- 
ginning of the century Bailey (1899, 1902) investigated the 
properties of RR Lyrae pulsators in the period-amplitude 
diagram which is now usually named after its inventor. 
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However, the investigation of pulsational properties in the 
Bailey diagram has never gone beyond an empirical and 
qualitative approach. For instance, we know that along 
this diagram the RR Lyrae variables of different Ooster- 
hoff types show different distributions, but until recently 
we have been unable to investigate the origin of such a 
behavior because we lack a detailed theoretical scenario 
connecting pulsational amplitudes with the structural pa- 
rameters of the pulsating stars. This is due to the fact 
that periods can be easily derived by a simple linear and 
adiabatic approach to the pulsational phenomenon, i. e. 
the canonical von Ritter relation, whereas the amplitudes 
are the result of detailed nonlinear investigations which re- 
quire a much more difficult theoretical and computational 
effort. 

To address this problem properly, during the last few years 
we have undertaken a theoretical project aimed at the 
evaluation of the pulsational amplitudes and the modal 
stability of RR Lyrae variables by computing a large and 
homogeneous set of nonlinear and time-dependent convec- 
tive envelope models. The primary goal of this approach is 
to provide, for the first time, a sound theoretical scenario 
of the dependence of all the pulsational observables on as- 
trophysical parameters such as stellar mass and chemical 
composition. Moreover, the pulsation properties of nonlin- 
ear RR Lyrae models, when evaluated at full amplitude, 
can supply useful suggestions concerning the physical pro- 
cesses which drive and/or quench the pulsation mecha- 
nism. 

Bono & Stellingwerf (1994, hereinafter referred to as BS) 
have already shown that a nonlinear approach is necessary 
for an accurate determination of periods. The outcome 
was similar for properly locating the instability regions 
into the HR diagram. Indeed, the nonlinear instability 
boundaries of fundamental (F) and first overtone (FO) 
pulsators have been already found in substantial agree- 
ment with the observed color distribution of RR Lyrae 
variables in selected galactic globulars (Bono, Caputo & 
Marconi 1995, hereinafter referred to as BCM). Moreover, 
the minimum fundamentalized periods of RR pulsators 
belonging both to Oosterhoff I (Ool) and Oosterhoff II 
(OoII) globular clusters have been explained by taking 
simultaneously into account the evolutionary history and 
the modal stability of variable stars in the instability strip 
(Bono et al. 1995). 

The results presented in this paper have been obtained by 
using the same theoretical framework adopted by Bono 
& Stellingwerf (1992) and BS in the first two papers of 
this series. Physical assumptions and numerical proce- 
dures have been already described in the quoted papers 
and in Stellingwerf (1975, 1982). The reader is referred 
to these papers for a detailed description. The main dif- 
ferences between the present investigation and the results 
obtained by BS are: 

1) the RR Lyrae models supplied by BS were computed 
by assuming a fixed mass value (M=0.65 M & ) whereas in 
the present paper three different mass values have been 
adopted (M=0.58, 0.65, 0.75 M@) which cover the whole 
expected range of HB stars. 

2) In order to compare the pulsation properties of RR 
Lyrae models with previous theoretical works, BS adopted 



a helium abundance by mass of Y=0.30 whereas all the 
present computations have been performed by adopting a 
helium content of Y=0.24. We chose this value since on the 
basis of both the R (or R') method (Buzzoni et al. 1983) 
and the determination of the ratio He/H in HII regions 
and planetary nebulae (Peimbert & Torres-Peimbert 1976; 
Peimbert 1995 and references therein) there is a general 
consensus on the fact that a suitable range for the helium 
abundance of population II stars belonging to GGCs is 
0.22 < Y < 0.25 (Pagel 1995 and references therein). 
Owing to the marginal dependence of the pulsational 
properties of metal poor stars on metal content (Christy 
1966), a fixed metal abundance, Z=0.0001, has been 
adopted. However, it is worth underlining that throughout 
the pulsation cycle the transformation of luminosity and 
temperature into the observational plane depends on the 
metal content, since both bolometric correction and color- 
temperature relation present a non negligible dependence 
on such a parameter. 

According to theoretical expectations (see, e.g., Castel- 
lani, Chieffi & Pulone 1989, hereinafter referred to as 
CCP) we investigated the pulsational behavior of clus- 
ter variables by assuming two different stellar masses for 
properly covering the range of pulsator masses expected in 
clusters with "solar scaled" metallicities Z~0.0001 (M = 
0.75M©) and Z~0.001 (M = 0.65M©). In order to account 
for the suggested overabundance of alpha-elements, which 
would produce an increase in the "effective metallicity" 
(Straniero & Chieffi 1991) and a further decrease in the 
pulsator masses, a set of models characterized by a smaller 
mass value, M = 0.58 M & , was also computed. According 
to CCP, this pulsator mass corresponds to a cluster effec- 
tive metallicity Z~0.006. This value can be obtained by 
increasing, in a solar scaled metallicity of Z=0.001, the 
amount of alpha-elements of a factor roughly equal to 9 
(see Salaris, Chieffi & Straniero 1991). 
In the following section we provide an atlas of selected 
theoretical light curves arranged according to the effec- 
tive temperature and the luminosity of each model. In §3 
nonlinear periods are compared to the periods obtained 
by using the van Albada and Baker (1971) relations. In 
order to accomplish a comparison with observational data 
from RR Lyrae stars belonging to GGCs, the results of 
pulsation theory will be connected with the canonical HB 
evolutionary scenario. 

Section 4 deals with the pulsational amplitude as a func- 
tion of effective temperature, luminosity and mass. These 
results are collected in §5 for providing a theoretical Bai- 
ley diagram which is examined in terms of evolutionary 
expectations. In §6 this theoretical scenario is compared 
with observational evidences from RR Lyrae pulsators in 
selected GGCs. A final discussion about some residual 
but relevant mismatches between theory and observations 
closes the paper. 

2. Light Curves 

For each given value of the stellar mass the modal stability 
of the envelope models has been explored at selected lumi- 
nosity levels (logL/L @ = 1.51, 1.61, 1.72, 1.81, 1.91) and 
by assuming a step of 100 K in the effective temperature. 
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Figure 1 shows the topology of the instability strip for 
M=0.58 Mq, as compared with similar boundaries but for 
the mass values M = 0.65 or O.75M0 already reported in 
Bono, Caputo & Marconi (1995). Table 1 gives luminosi- 
ties and effective temperatures along the newly computed 
boundaries. One finds that the topology of the instability 
strip for the mass value M = O.58M0 discloses the same 
overall features of the instability strip of RR Lyrae models 
characterized by larger mass values. From the same figure 
it can be seen quite clearly that the width of the region 
where only the FO presents a stable limit cycle increases 
as the stellar mass increases. 

Figure 1 also shows the location of some FO models for 
M = O.58M0 which present a stable limit cycle at effective 
temperatures lower than the FO red edge. These models 
define a FO stability isle in a region of the instability strip 
where only F pulsators should be characterized by a stable 
limit cycle. 

A similar but specular finding was obtained by BS. Indeed, 
they found some F models showing a stable limit cycle in 
the region of the instability strip where only FO pulsators 
might be present. The approach to limit cycle stability 
of these models due to the peculiarities of the dynamical 
and convective structure will be discussed in a forthcoming 
paper (Bono, Marconi & Stellingwerf 1996). 



Fig. 1. The location into the HR diagram of the instability 
strip for M = 0.58 M® is compared with previous results for 
larger values of stellar mass and fixed chemical composition. 
At the lower luminosities, and for each assumed mass value, 
moving from higher to lower effective temperatures the differ- 
ent curves show: the FO blue edge (FO-BE), the F blue edge 
(F-BE), the FO red edge (FO-RE) and the F red edge (F-RE). 
Asterisks mark the location of some models for M = 0.58 M® 
which are characterized by a stable first overtone limit cycle 
in the region of the instability strip where only fundamental 
pulsators should be present. See text for further details. 



Table 1. Fundamental and first overtone instability bound- 
aries for M = O.58M . 



logi/L FO-BE F-BE FO-RE F-RE 





(°K) 


(°K) 


(°K) 


(°K) 


1.86 


6850 




6150 




1.81 


6950 


7050 


6750 


5750 


1.77 


7050 


7050 


6750 




1.72 


7150 




6650 




1.61 


7250 


7150 


6750 


5950 
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In order to achieve a good accuracy throughout all phases 
of the pulsation cycle, a variable time step has been 
adopted: the number of time steps necessary for covering 
a period ranges roughly from 400 to 600. Before being able 
to provide any sound conclusion concerning the modal sta- 
bility and the dependence of the light and velocity curve 
morphologies on physical parameters and chemical com- 
position, the pulsation characteristics have to approach a 
limiting amplitude. To accomplish this goal the local con- 
servation equations and the convective transport equation 
for each case were integrated in time until the initial veloc- 
ity perturbation settles down and the pulsation behavior 
(amplitudes, periods) approaches a limit cycle stability. 
The time interval necessary for the radial motions to ap- 
proach their asymptotic amplitudes ranges from one thou- 
sand to ten thousand pulsational cycles. As a general 
rule, we assumed the direct time integration to be roughly 
equivalent to the inverse of the linear growth rate. There- 
fore, the nonlinear unstable models of the present survey 
are all characterized, over two consecutive periods, by a 
periodic similarity of the order of or lower than 10~( 4 ~ 5 ). 
The range of luminosities covered by the different series 
of models was chosen so as to largely cover the theoreti- 
cal prescriptions from HB evolutionary models. For vari- 
able stars these models foresee a minimum Zero Age Hor- 
izontal Branch (ZAHB) luminosity log Lzahb — 1-70 for 
Z=0.0001 and log L Z ahb ~ 1-65 for Z=0.001 respectively. 
At the same time, the canonical HB evolutionary scenario 
predicts the occurrence of variable stars which are evolving 
within the instability strip at moderately larger luminosi- 
ties. Figure 2 presents an atlas of theoretical light curves 
of F pulsators for two consecutive periods. These models 
were computed by assuming a mass M=0.65 M @ . The light 
curves of FO pulsators referred to the same sequences of 
models are reported in Figure 3. Figures 4 and 5 show the 
same quantities of Figure 2 and 3 but are referred to the 
sequences in which a mass value M=0.75 M @ was adopted. 
Both light and velocity curves for the models discussed in 
this paper are available upon request to the authors. 



Fig. 2. Theoretical light curves of the four sequences of funda- 
mental pulsators for the stellar mass M=0.65 M©. Each plot 
shows the bolometric amplitude for two consecutive periods 
and the effective temperature of the model. The luminosity 
level of the four sequences are also reported. 



Fig. 3. Same as Fig. 2, but the light curves are referred to the 
four sequences of first overtone pulsators. 
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Fig. 4. Theoretical light curves of the four sequences of funda- 
mental pulsators for the stellar mass M=0.75 M©. Each plot 
shows the bolometric amplitude for two consecutive periods 
and the effective temperature of the model. The luminosity 
level of the four sequences are also reported. 



Fig. 5. Same as Fig. 4, but the light curves are referred to the 
four sequences of first overtone pulsators. 

From the calculated light curves we find that the overall 
morphology, as well as the occurrence of secondary fea- 
tures like "Bumps" or "Dips" over the pulsation cycle 
closely follows the theoretical scenario suggested by BS 
which is based on models computed by adopting a differ- 
ent helium abundance (Y=0.30). The overall agreement 
concerning the occurrence of secondary features both in 
theoretical and observational light curves has been already 
discussed by Walker (1994) by presenting the photometry 
of RR Lyrae stars in the galactic globular M68. 

3. Nonlinear periods 

Figures 6 and 7 show the distribution of periods for all 
computed models with mass M=0.65 and 0.75 M & which 
reach a stable limit cycle either in the F or in the FO mode. 
In the same figures we report the linear predictions given 
by van Albada & Baker (1971) according to the relations: 

logP F = 11.497 + 0.84 log L- 0.68 log M - 3.48 log T e 

log(-^-) = 0.43 + 0.14 log L — 0.032 log M — 0.009 log T e 

Pfo 

where both mass and luminosity are in solar units, whereas 
the periods are in days. 



Fig. 6. Fundamental periods versus effective temperatures for 
two different stellar masses M/Mq =0.65 (open circles) and 
M/Mq =0.65 (solid circles). Solid lines refer to the nonlinear 
results of the present computations, whereas dotted lines show 
the linear periods evaluated by using the relations provided 
by van Albada & Baker (1971). The luminosity level of each 
sequence is also reported at the right of the curves. 

As expected, nonlinear periods are generally smaller than 
linear periods, since convective motions smooth the den- 
sity profile in coincidence of the hydrogen ionization region 



Fig. 7. Same as Fig. 6, but for first overtone pulsators. 

where the density inversion takes place, and in the mean 
time produce a non negligible increase of the density in the 
region where the adiabatic exponent -Ti- reaches its min- 
imum value. The aftermaths of this effect together with 
the nonlinear effects on Ti cause the decrease of the non- 
linear periods when compared to linear ones. The previous 
leading-term considerations also provide an explanation of 
the physical reason why this discrepancy tends to increase 
moving from the blue to the red side of the instability strip 
(see Figure 11 and 12 in BS). 

The new models show that the differences are larger for 
FO than for F pulsators, whereas BS found that the dis- 
crepancy between linear and nonlinear periods was larger 
for F pulsators. Since the present nonlinear models differ 
from the models adopted in BS only in the helium abun- 
dance, this finding strongly suggests that the FO periods 
present a non negligible dependence on this parameter. 
Linear interpolation among the data plotted in both fig- 
ures gives the new relations: 

logP F = 11.627 + 0.823 log L- 0.582 log M- 3.506 log T e 
log Ppo = 10.789 + 0.800 log L — 0.594 log M — 3.309 log T e 

with a correlation coefficient r=1.00. The differences be- 
tween the previous relations and the van Albada & Baker 
relations -though not negligible- are not expected to pro- 
duce a substantial change in the available pulsational sce- 
nario, supporting once more the pioneering work by van 
Albada & Baker. However, previous results give at least 
a firm warning against the long-standing attempt to de- 
rive the mass of double-mode pulsators by using the ratio 
between F and FO periods {Pf I Pfo)- Indeed, this quan- 
tity critically depends on both stellar parameters and on 
even minor details of the physical and numerical ingre- 
dients adopted in the pulsational codes. The occurrence 
and the properties of RR Lyrae models which show stable 
double-mode limit cycles will be discussed in a forthcom- 
ing paper. 

The previous relations, which provide the nonlinear peri- 
ods for the first two pulsating modes, can be combined 
with the instability strip topologies given by BCM to de- 
rive the behavior of the minimum period allowed in a clus- 
ter, i. e. the minimum (fundamentalized) period of the FO 
pulsator at the blue edge of its instability region: 

log Pfo-be = -2.424 - 0. 714 log M + 1.046 log L 

which appears in good agreement with the preliminary 
evaluation obtained by Bono et al. (1995). As discussed in 
the quoted paper the minimum periods provided by the 
previous relation reproduce quite satisfactorily the mini- 
mum periods of variable stars belonging to GGCs, thus 
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supporting once more the canonical evolutionary scenario 
against the suggested occurrence of an anomalous period 
shift. 

As a further step in our investigation, let us take now into 
account the canonical evolutionary constraints given by 
CCP. Figure 8 shows the expected location of HB evolu- 
tionary tracks for selected values of the cluster metallicity. 
This parameter was chosen such that it should become 
representative of well known clusters which populate the 
galactic halo and present a large number of RR Lyrae 
variables, namely Z=0.0001 (M15), Z=0.0004 (M3) and 
Z=0.001 (M5). 

In the same figure we plot the instability strip for the two 
labeled assumptions about the mass of the star, giving, 
for each assumed metallicity the maximum mass value 
(Mi£ ax ) allowed by evolutionary models for populating 
the HB of a cluster characterized by an age of 15 Gyr and 
the expected minimum period of ab type pulsators (see 
below). It appears that the instability strip in M15 should 
be populated by stars with mass values included between 
0.75 and 0.80 M @ . More massive stars would be forbidden 
for the adopted cluster age, thus giving a straightforward 
explanation of the blue HB shown by this cluster. 
On the contrary, the instability regions of M3 and M5 
should be populated by stars with mass values which range 
from 0.65 to 0.70 M & respectively. These values imply now 
an amount of mass loss of the order of 0.10 M @ or 0.18 
M & for M3 and M5 respectively. 



This interesting observational behavior has stimulated an 
animated debate during the past few years. According to 
Sandage (1958, 1981, 1982, 1990) the difference in periods 
should mainly be ascribed to a difference in luminosity: 
OoII pulsators are more luminous than Ool ones. 
In order to explain the F mean period distribution, van Al- 
bada & Baker (1973) suggested an alternative hypothesis: 
a difference in the maximum temperature of F pulsators 
produced by a hysteresis mechanism. In this context a 
variable star located in the OR region will pulsate in the 
F mode if it is evolving from lower to higher effective tem- 
peratures, whereas it will pulsate in the FO if the star is 
evolving in the opposite direction. As a consequence of this 
mechanism, we expect that the transition between RR a j 
and RR C occurs either at the F blue edge or at the FO red 
edge. 

In the already quoted Figure 8, the point where the tran- 
sition is expected to take place according to the hystere- 
sis hypothesis is marked with an open circle. Accordingly 
one expects that the OR region is populated by FO pul- 
sators in OoII clusters, and by F pulsators in the Ool case. 
According to such an hypothesis, the data listed in Ta- 
ble 2 disclose that the theoretical results concerning P™ n 
overlap surprisingly well with the observational data when 
both evolutionary prescriptions and pulsational results are 
taken into account. 

If one adds the evidence given in BCM for a HR distri- 
bution in agreement with the hysteresis prescriptions, we 
can conclude that the canonical evolutionary scenario ap- 
pears able to account for several relevant features of RR 
Lyrae pulsators, without the need of invoking "ad hoc" 
modifications of the current evolutionary knowledge. 



Fig. 8. The ZAHB location into the HR diagram for three 
different values of the cluster metallicity. In each plot the solid 
lines show the evolutionary paths of selected helium burning 
models for the labeled values of the masses. Dashed lines and 
dotted lines report the instability strip for two different mass 
values M=0.75 M and M=0.65 M respectively. For each 
given metallicity we also report the name of the typical clus- 
ter, the observed minimum period of fundamental pulsators 
in that cluster (logP™"), and the maximum mass allowed for 
populating the HB (M^ aa: ). See text for further explanations. 

As a most relevant point, we find that this evolutionary 
picture can be easily connected with the observed mini- 
mum period of RR a j pulsators, i.e. with the observational 
evidence at the basis of the Oosterhoff dichotomy. The 
first suggestion about the occurrence of the Oosterhoff di- 
chotomy was produced by the observational evidence that 
the mean periods of F pulsators in GGCs, < P a b >, were 
grouped around two distinct values, namely 0.65 d (OoII) 
and 0.55 d (Ool), with alack of clusters with < P a b >=0.6 
d. As it was early recognized, such an occurrence is due to 
the fact that when compared to OoII clusters, Ool clus- 
ters allow the occurrence of F pulsators characterized by 
shorter periods. Indeed, the minimum F period changes 
from log P% in ~ -0.25 in OoII clusters to log P% in ~ - 
0.35 in Ool clusters. 



Table 2. Theoretically predicted minima of fundamental pe- 
riods as a function of cluster metallicity. 



z 


M 


logi 


logT e 


logP a f" 


Transition 


0.0001 
0.0004 
0.0010 


0.80 
0.70 
0.65 


1.761 
1.710 
1.659 


3.820 
3.844 
3.843 


-0.26 
-0.35 
-0.37 


FO-RE 

F-BE 

F-BE 



As a final point, let us note that the data in Table 2 give 
a spontaneous indication for the shift in periods within 
the Ool type clusters (Castellani & Quarta 1987), with F 
pulsators in M5 reaching shorter periods with respect to 
the same type of pulsators in M3. 

4. Pulsational amplitudes 

A first attempt to connect pulsational amplitudes with 
stellar parameters has been recently presented by Bro- 
cato, Castellani & Ripepi (1996) who discussed new data 
for RR Lyrae stars in the globular M5. However, their 
investigation was necessarily based on the available pulsa- 
tional computations given by BS for a fixed but unrealistic 
amount of atmospheric helium (Y=0.30). Similar pioneer- 
ing results are now being compared with amplitudes com- 
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puted under the much more realistic assumption Y=0.24 
(Pagel 1995, Peimbert 1995). 

Figure 9 shows the bolometric amplitude as function of 
the effective temperature for a star of 0.65 M & at selected 
stellar luminosity values and for the two labeled assump- 
tions of atmospheric helium Y=0.24 or Y=0.30. Note that 
the figure reports only the amplitudes of models which 
present a stable limit cycle in the F and/or in the FO. 
One finds that the relation between period and tempera- 
ture obtained by assuming a helium content Y = 0.24 does 
not differ considerably from that obtained by assuming a 
larger helium abundance (Y = 0.30). The small differences 
in the range of temperatures taken into account appear as 
the direct consequence of the different location in the HR 
diagram of the corresponding instability regions. 



Fig. 9. Bolometric amplitude versus effective temperature for 
both fundamental and first overtone variables. The amplitudes 
plotted in this figure are referred to models computed by as- 
suming a fixed stellar mass M=0.65 and two different 
helium abundances Y=0.24 (solid lines) and Y=0.30 (dashed 
lines). The sequences of models characterized by different lu- 
minosity levels are marked with different symbols. 

A relevant feature for Y=0.24 is that the correlation be- 
tween the amplitudes of F pulsators and the effective tem- 
perature becomes more and more independent of stellar 
luminosity in comparison with the models which present 
a higher helium abundance (Y=0.30). This is of course 
a crucial point, since it leads to a direct correlation be- 
tween F-amplitudes and effective temperatures, as earlier 
suggested by Sandage (1981). Moreover, one finds that 
the distributions of FO amplitudes versus effective tem- 
perature keeps showing the characteristic "bell" shape, 
with the maximum amplitude sensitively dependent on 
the luminosity level of the pulsators. These amplitudes 
moving from Y=0.30 to 0.24 do not show systematic dif- 
ferences, since at higher luminosities the two curves are 
almost identical. At lower luminosities the models with 
Y=0.30 present pulsational amplitudes which are initially 
larger and then smaller in comparison with the models 
characterized by Y=0.24. 

Before attempting to match pulsational and evolutionary 
prescriptions, it is necessary to evaluate the effects of vary- 
ing the stellar mass in the above quoted scenario. This is 
shown in Figure 10 in which we report pulsational ampli- 
tudes against effective temperatures but for fixed helium 
content (Y=0.24) and for two different assumptions on 
stellar mass M=0.65 and 0.75 M@. 

Figure 10 shows that a mass increase has little effect 
on the amplitudes of F-pulsators, whose independence of 
the luminosity level appears even strengthened for higher 
masses. On the contrary, at lower luminosities we find that 
an increase of the stellar mass implies an increase in FO 
amplitudes. As a point that will be relevant further on, 



Fig. 10. Bolometric amplitude versus effective temperature 
for both fundamental and first overtone variables. The ampli- 
tudes plotted in this figure are referred to models computed 
by assuming a fixed helium content (Y=0.24) and two differ- 
ent stellar masses M=0.65 M (solid lines), and M=0.75 M 
(dashed lines). The sequences of models characterized by dif- 
ferent luminosity levels are marked with different symbols. 



Figure 10 supports the suggestion of BS that the ampli- 
tude of FO pulsators along the decreasing branch leading 
to the red limit is largely independent of the stellar lumi- 
nosity. We further find that, for a given temperature, this 
amplitude appears moderately dependent on the pulsator 
mass, increasing when the mass increases. We can find 
a reason for such a behavior by recognizing that the de- 
creasing branch appears populated by FO pulsators in the 
OR-region, and by assuming that the pulsational ampli- 
tude only depends on the distance in temperature from the 
red edge. According to such a picture, FO amplitudes are 
expected to be largely independent of the luminosity for 
a given mass, while increasing the mass at any given tem- 
perature increases the amplitude because the FO bound- 
ary shifts toward cooler temperatures, thus increasing the 
distance of the given temperature from the FO-RE. 

Concerning the predicted behavior of stars in GGCs, one 
may notice that the amplitude distributions shown in Fig- 
ure 10 give two different periods in the OR-region, one 
along the F and the other along the FO sequence. If the 
suggestion advanced in the previous section about the effi- 
ciency of the hysteresis mechanism is correct, one expects 
OR-regions populated according to the Oosterhoff type, 
either by FO (OoII) or by F pulsators (Ool). In terms of 
the data shown in Figure 10 one expects in OoII type clus- 
ters the occurrence of a well developed decreasing branch 
of FO pulsators together with a sequence of F pulsators 
depopulated in its upper portion. On the contrary, in Ool 
type clusters one expects a sequence of F pulsators popu- 
lated up to the largest amplitude, with negligible evidence 
for the decreasing branch of first overtones. 

The discussed scenario can be safely extended to the case 
for M=0.58 M@, as illustrated by data in Figure 11. In- 
deed, both the FO amplitudes and the width in temper- 
ature of the FO region regularly decrease when passing 
from M=0.75 M@ to 0.65 and 0.58 M @ . Moreover, one 
finds that F pulsators arrange even better along the pre- 
vious amplitude-temperature relation. As a result, if pul- 
sational prescriptions given through figures 10 and 11 are 
taken at their face values, the observed amplitudes of F 
pulsators should allow the determination of stellar effec- 
tive temperatures with an error of about ±100 K, if not 
better, independently of the stellar mass and luminosity. 
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Fig. 11. Pulsational amplitudes versus effective tempera- 
ture for three different assumptions on stellar mass. First 
overtone amplitudes refer to the common value of luminos- 
ity logL/Lo = 1.61. Fundamental amplitudes are given for 
logL/L @ = 1.61, 1.81 (M=0.58 M@) and logL/L @ = 1.72, 
1.81. 1.91 (M=0.65 M@). 

5. The theoretical Bailey period-amplitude dia- 
gram 

Theoretical periods and theoretical amplitudes, as ob- 
tained for each given assumption about stellar mass, 
can be finally collected to produce a theoretical period- 
amplitude diagram, thus allowing to investigate the con- 
nection of these two pulsational parameters with the evo- 
lutionary parameters. Figure 12 shows the distributions 
into this theoretical Bailey diagram when the stellar tem- 
perature is varied across the instability region, for selected 
values of the stellar luminosity and for the two choices 
M=0.65 M@ (solid lines) and M=0.75 M@ (dashed lines), 
keeping everywhere Y = 0.24. The pulsational amplitudes 
plotted in this figure are referred to four different lumi- 
nosity levels and cover the overall instability region. 



Fig. 12. Theoretical Bailey diagram, bolometric amplitude 
versus the logarithm of the period for a fixed helium content 
(Y=0.24) and two different mass values M=0.65 M (solid 
lines), M=0.75 (dashed lines). Symbols concerning the lu- 
minosity levels are the same as in Fig. 10. 

Bearing in mind the distribution reported in Figures 10 
and 11, the topology of data plotted in Figure 12 can 
be easily understood in terms of the relation which con- 
nects the periods with effective temperatures, luminosities 
and masses. By using these relations, for a given value 
of the mass and for a fixed luminosity level, the period 
scales according to log P oc 3.5 (or 3.3) log Te and the 
distributions shown in Figure 10 can be therefore trans- 
lated into the distributions shown in Figure 12. Moreover, 
by increasing the luminosity level the periods increase by 
about 8 log P ~ 0.8 log L. Such an occurrence removes the 
"degeneracy" of F amplitudes with stellar luminosity and 
produces the period shift observed in Figure 12 for a fixed 
value of the mass. According to these simple arguments, 
we finally find that the amplitude distribution, by increas- 
ing the mass, moves toward shorter periods, thus giving 
full account of the Bailey diagram morphology disclosed 
in Figure 12. 

Adopting such a pulsational scenario, it is obviously in- 
teresting to match pulsational constraints with the pre- 
dictions of the evolutionary theory. For a quick look on 
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this new scenario let us first discuss the distribution of 
HB stars in the Bailey diagram by assuming that all stars 
are just in their ZAHB location. Under this assumption, 
for each given value of the metallicity we derived the set 
of stellar masses populating the instability strip by in- 
terpolating the evolutionary tracks given by CCP, each 
mass being characterized by proper values of both effec- 
tive temperature and luminosity. The regular behavior of 
F pulsators allows a straightforward linear interpolation of 
the dependence of the amplitude on stellar mass, temper- 
ature and luminosity. On this basis evolutionary data can 
be easily transformed into the Bailey diagram for F pul- 
sators as given in Figure 13 for the three different choices 
of metal abundance Z=0.0001, Z=0.0004 and Z=0.001. 
The procedure previously outlined supplies a relevant re- 
sult that can be summarized as follow: 
the difference of the evolutionary structure of the pulsators 
predicted on the basis of different assumptions about the 
stellar metallicity plays a minor role in the distribution of 
the pulsators in the Bailey diagram. 

As suggested by Brocato, Castellani & Ripepi (1996) on 
a pure observational basis, this is due to the fact that 
moving from metal poor OoII to metal rich Ool pulsators 
the decrease in mass is largely balanced by the decrease 
in luminosity as can be seen in Figure 8. 
Unfortunately, the distribution of FO pulsators appears 
much less predictable. Due to the rather intricate ampli- 
tude behaviors disclosed in Figures 10 and 12, interpo- 
lation over the whole range of theoretical data appears 
not every useful, until a much finer set of nonlinear lumi- 
nosity amplitudes becomes available. However, the regular 
behavior of the decreasing branch allows again a linear in- 
terpolation of the dependence of such a feature on stellar 
parameters. On this basis the distribution of FO pulsators 
could be predicted. 

We eventually underline that the distribution of F pul- 
sators in the Bailey diagram, even when taking into 
account the off ZAHB evolution, would be affected 
only marginally, since these variables are characterized 
by a well defined amplitude-temperature relation which 
presents a negligible dependence on both stellar mass and 
luminosity level (see Figure 9, 10 and 11). The conclusions 
previously drawn for F pulsators cannot be extended to 
FO pulsators due to the strong dependency of their pul- 
sational amplitudes on luminosity level. As a consequence 
the distribution of these variables in the Bailey diagram, 
with the exception of pulsators located along the decreas- 
ing branch, might be partially changed by the off ZAHB 
evolution. 



Fig. 13. The distribution in the theoretical Bailey diagram 
of ZAHB fundamental and first overtone pulsators as expected 
according to evolutionary prescriptions for the labeled assump- 
tions about the stellar metallicity. The distribution of the de- 
creasing branch of first overtone pulsators is also sketched. See 
text for further explanations. 
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6. Comparison with observations 

To approach the problem of comparing the present theo- 
retical scenario with observational data, in Figure 14 we 
show the observational Bailey diagram given by Sandage 
(1990, and references therein) for pulsators belonging to 
the two well studied clusters M15 (OoII) and M3 (Ool). 
A quick look to the distributions plotted in this Figure 
reveals a series of relevant similarities with theoretical re- 
sults that are worth listing: 

i) F pulsators in both Ool and OoII clusters arrange along 
a rather continuous sequence. 

ii) The OoII prototype M15 lacks F pulsators in the up- 
per portion of the amplitude distribution, whereas a large 
amount of M3 (Ool prototype) RR a b variables populate 
this region. 

iii) The lack of large amplitude F pulsators in M15 is ac- 
companied by the clear evidence for a decreasing branch 
of FO pulsators. On the contrary, M3 shows a much less 
clear evidence of such a decreasing branch, if any. 



Fig. 14. The observational Bailey diagram, blue amplitude 
versus the logarithm of the period, for RR Lyrae variables be- 
longing to galactic globulars M3 (solid circles) and M15 (open 
circles) as given by Sandage (1990). 



According to previous findings, we can conclude that the 
proposed theoretical scenario offers an interesting clue to 
an evolutionary interpretation of the observed features 
shown by the Bailey diagram supporting the occurrence 
of F (or FO) pulsators in Ool (OoII) clusters. However, 
a more detailed comparison reveals several quantitative 
mismatches which will be discussed in detail. 
To allow such a quantitative comparison with actual clus- 
ter pulsators, theoretical light curves have been trans- 
formed through Kurucz's (1992) evaluation of bolometric 
corrections and color-temperature relations to derive am- 
plitudes in the B photometric band, which is the most 
adopted observational parameter for RR Lyrae ampli- 
tudes. As it is well known, when the luminosity of the pul- 
sator decreases the star becomes cooler and the bolometric 
correction increases. As a consequence, B amplitudes for 
RR Lyrae pulsators are expected to be larger than the cor- 
responding bolometric amplitudes, an occurrence which 
of course justifies the widespread use of such an observa- 
tional parameter. Panel a and b of Figure 15 show the pre- 
dicted amplitude-temperature and amplitude-period rela- 
tions transformed in terms of B magnitudes. 
On the basis of this new pulsational topology the analysis 
given in Figure 13 can now be repeated for providing the- 
oretical predictions to be compared with observations. We 
report in Table 3 for the metal content Z=0.0001 the ex- 
pected pulsational properties of F and FO pulsators along 
the ZAHB. 



Fig. 15. The theoretical (a) amplitude-temperature and (b) 
amplitude- period relations for RR Lyrae variables as given in 
terms of blue amplitudes A(B). 



Table 3. Predicted ZAHB fundamental and first overtone pul- 
sators for Z=0.0001 (mass and luminosity are in solar units). 
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If we take into account only the decreasing branch of the 
FO amplitudes displayed in the left panel of Figure 15, for 
not too large luminosities, we derive the following relation: 



A{B) = 37.73 log T e + 1.93 log M - 143.29 

On the basis of such a relation we can attempt a predic- 
tion for variable stars located in the OR region, i.e. for 
stars already evaluated as F pulsators at the top of the 
F-branch. Data for these FO pulsators are also reported 
in Table 3. 

The comparison with observational data for pulsators be- 
longing to M15 is shown in Figure 16. In the same fig- 
ure are also plotted the evolutionary paths of F pulsators 
evolving off their ZAHB loci. Comparison with data for 
M15 discloses a surprising agreement. However, the same 
figure shows that theory foresees the occurrence of large 
amplitudes (A(B)~1.1 mag) for FO pulsators which are 
not present in the cluster. 



Fig. 16. Observational data for M15 pulsators are compared 
with the theoretical distribution for Z=0.0001 ZAHB pulsators. 
In the same figure are also shown the evolutionary paths for 
the labeled values of stellar masses. 
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The procedure has been repeated adopting Z=0.0004, as- 
sumed as a suitable choice for Ool clusters like M3. Data 
for both F and FO pulsators are reported in Table 4, 
whereas in Figure 17 theoretical predictions are compared 
with observational data for the quoted clusters. The agree- 
ment appears now somehow less satisfactory, since some 
F pulsators tend to have lower periods than predicted by 
the theory. If this occurrence is not an artifact of forcing 
the results for the lower luminosity levels into a common 
linear interpolation, we find no obvious origin for such a 
behavior. 



Table 4. Predicted ZAHB fundamental and first overtone pul- 
sators for Z=0.0004. 



about 300 K larger than predicted by the theory. Note 
that data in Figure 18 assume a reddening value E(B- 
V)=0.07 for M15 (Zinn 1985). By adopting the reddening 
value suggested by Sandage (1990), i. e. E(B-V)=0.11, 
the discrepancy for this cluster further increases with 200 
K. The reader is referred to Brocato, Castellani & Ripepi 
(1996) and to Silbermann & Smith (1995) for a detailed 
discussion concerning the M15 reddening evaluations. 
As a conclusion, the theoretical pulsational scenario we 
are dealing with, although it presents relevant similarities 
with the actual behavior of cluster pulsators, is far from 
reaching a tight correspondence with the available obser- 
vational data. 
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As for F pulsators once again we find that theoretical am- 
plitudes for FO pulsators are larger than observed values. 
However, we look at the comparison in Figure 17 as an ev- 
idence of the lack of the decreasing branch, disregarding 
the small amplitude variations at the larger periods as an 
evidence for the present location of the decreasing branch 
phase. 



Fig. 17. Observational data for M3 pulsators are compared 
with the theoretical distribution for Z=0.0004 ZAHB pulsators. 
In the same figure are also shown the evolutionary paths for 
the labeled values of stellar masses. 



As a final point, in Figure 18 are compared theoretical 
and observational data concerning the dependence of blue 
amplitudes on stellar effective temperatures. According to 
Brocato, Castellani & Ripepi (1996) the pulsator mean 
colors have been first corrected for reddening and then 
transformed into the true colors of the static models, us- 
ing the procedure given by Bono, Caputo & Stellingwerf 
(1995). Temperatures have been finally derived adopting 
the color-temperature relation as given by Kurucz(1992). 
We find a further disagreement between observations and 
theory. Indeed observational amplitudes appear smaller 
for each given value of the effective temperature or, con- 
versely, for each given amplitude the temperature appears 



Fig. 18. Theoretical expectations concerning the B ampli- 
tude-temperature relations for Z=0.0001 or Z=0.0004 are com- 
d with observational data for the cluster M15 and M3, re- 
spectively. Equilibrium temperatures are from Bono, Caputo 
& Stellingwerf (1995). 



One could take the comparison in Figure 18 as an evidence 
that our theoretical approach could be affected by a sys- 
tematic overestimate of amplitudes. As a matter of fact, 
the difference in temperature of about 300 K for a given 
amplitude appears much larger than the usual estimates of 
errors in the color-temperature relation of pulsating stars. 
If this is the case, the surprising agreement we found for 
the Bailey diagram of metal poor clusters should be an 
artifact of balancing errors. 

However, in order to properly set the present observational 
scenario concerning the systematic errors in the RR Lyrae 
effective temperature scale, it is worth mentioning that, 
by adopting the V-R color-temperature relation provided 
by Vanderberg & Bell (1985) instead of Kurucz's relation, 
Silbermann & Smith (1995) obtained a systematic shift 
toward cooler temperatures. Most interesting, they also 
found that by adopting the temperature calibration de- 
rived by Longmore et al. (1990) on the basis of V-K colors, 
the RR Lyrae temperatures are in average 300 K cooler 
than the temperature obtained by using the V-R colors. 
Here we can only present this contradictory scenario as 
something to be explored before reaching the final goal 
of reading the location of RR Lyrae pulsators in the Bai- 
ley diagram in terms of the evolutionary status of cluster 
stars. 



7. Discussion and conclusions 

Before closing the argument, let us discuss if and how 
the theoretical scenario could be varied in the attempt to 
match the theoretical results with observations. The theo- 
retical framework outlined in the present paper is based on 
numerical experiments supplied by three different fields: 
1) stellar evolution, 2) stellar pulsation, 3) stellar atmo- 
spheres. Based on proper physical assumptions, all of them 
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simulate the physical processes which take place in differ- 
ent regions of the stellar structure. Unfortunately, any ob- 
servational test involves in the same context the results of 
the quoted scenario so that it appears difficult to point out 
where and why the theoretical outcomes could be affected 
by systematic errors. In order to give some more light on 
this scenario, let us briefly discuss in the following the 
three above quoted theoretical ingredients. 

1) Several and thorough papers have been recently de- 
voted to the evolutionary properties of HB stars and to 
their dependence on astrophysical parameters and phys- 
ical ingredients (CCP; Dorman, Rood & O'Connell 1994 
and references therein) and therefore they will not be dis- 
cussed here. 

2) The theoretical treatment of radial stellar pulsations is 
based on two main ingredients: the physical inputs nec- 
essary to produce an envelope model and the set of con- 
vective hydrodynamic equations adopted for following the 
nonlinear coupling between dynamical and convective mo- 
tions. 

Among the physical inputs, both radiative opacities and 
equations of state are worth being checked, since the ra- 
dial pulsational instability is driven by the destabilizing 
effects operated by kappa and gamma mechanisms. Nu- 
merical tests performed to verify the dependence of the 
nonlinear pulsation characteristics on the new OP (Seaton 
et al. 1994) and OPAL (Roger & Iglesias 1992) opaci- 
ties show almost no difference in comparison with models 
based on the old Los Alamos opacities. A similar result 
was obtained by Buchler & Buchler (1994) in their inves- 
tigation on BL Herculis variables. Therefore this ingredi- 
ent, at least for low metallicity stars, has no effect on the 
amplitude problem. 

Concerning the equation of state we did not perform any 
specific calculation since recent nonlinear results on radia- 
tive models of bump Cepheids provided by Kanbur (1992) 
show only negligible physical differences on the adopted 
equations of state. However, we plan to test the new Liv- 
ermore equation of state as soon as it becomes available 
for metal contents suitable to metal poor stars. 
Although the assumptions we adopted to describe the tur- 
bulent field in variable stars are physically plausible, the 
convective transport equation has been derived by using 
three different free parameters (for a detailed discussion 
see Stellingwerf 1982 and BS). A more formal and rig- 
orous approach for deriving the set of radial pulsation 
equations was adopted by Kuhfufi (1986) and by Gehmeyr 
(1992 and references therein). However, only few full am- 
plitude nonlinear models have been computed by adopting 
the latter formalism and therefore a trustworthy compar- 
ison between hydrocodes which involve different physical 
and numerical approximations cannot be provided. 
The self-consistent inclusion of nonlocal, nonlinear and 
time-dependent effects may help to improve long-standing 
pulsation questions (BS; BCM; Bono, Caputo & Stelling- 
werf 1995). In particular, the nonlinear convective models 
are characterized by pulsation amplitudes that, although 
still too large, are roughly half magnitude smaller than 
the amplitudes of radiative models. 

3) The transformation of the light curves into the observa- 
tional plane is one of the most thorny problems of radial 



stellar pulsations. At present two different routes can be 
followed to accomplish this transformation. The former 
involves the calculation of hydrodynamic pulsating atmo- 
sphere models which take into account a multifrequency 
radiation transport in the outermost optically thin layers. 
This approach not only provides useful insights on the 
dynamical properties of the surface regions but can also 
predict broad-band colors. Until now only few investiga- 
tions have been devoted to soundly solve this problem for 
RR Lyrae variables (Keller & Mutschlecner 1970; Bendt 
& Davis 1971). 

Davis & Cox (1980) following this route computed few 
models for evaluating the RR Lyrae mean colors but no 
theoretical study has been undertaken so far for estimating 
the RR Lyrae light curves in specific photometric bands. 
The latter route relies on the use of bolometric corrections 
and color-temperature relations provided by static atmo- 
sphere models. This approach is based on two weak physi- 
cal assumptions: the stellar effective temperature over the 
full cycle is derived by assuming that the surface zone is al- 
ways in radiative equilibrium, whereas the surface gravity 
is evaluated by assuming both radiative and hydrostatic 
equilibria (Bono, Caputo & Stellingwerf 1994). Although 
the quoted approximations were proved to be suitable for 
evaluating the RR Lyrae mean colors (Bono, Caputo & 
Stellingwerf 1995) and for reproducing the ionization and 
excitation equilibrium conditions of metal lines close to 
the phase of minimum light (Clementini et al. 1995), the 
evaluation of pulsation amplitudes may require a more 
cautious analysis. 

We conclude that the present theoretical scenario proves 
to be successful to shed light on many pulsation features 
of cluster variables and, in particular, on long-standing 
questions such as the transition between FO and F pul- 
sators and the origin of the Oosterhoff dichotomy. How- 
ever, predictions concerning the pulsation amplitudes are 
far from reaching the required compatibility with ob- 
servational data. The possible origins of such a discrep- 
ancy were discussed, but firm conclusions have not been 
reached. 

We can only draw attention on such a situation, wait- 
ing for further improvements either of the theory of radial 
pulsation or of the theory of stellar atmospheres and/or of 
stellar evolution which may succeed in reconciling theory 
with observation. This will reach the relevant goal of using 
the pulsational parameters of RR Lyrae variables to gain 
firm information on their evolutionary status and, in par- 
ticular, on their masses, temperatures and luminosities. 
The use of these pulsators will allow us to solve many open 
questions concerning the evolution of our Galaxy and, in 
particular, they might become robust standard candles for 
mapping the distances inside the Galaxy and among the 
galaxies belonging to the Local Group. 
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